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ABSTRACT 

The  use  of  superconductors  in  the  field  windings  of  large 
synchronous  rotating  electrical  machines  virtually  eliminates  the 
need  for  magnetic  iron  core  material  and  hence  reduces  the  size  and 
weight  requirements  for  a  given  power  level.  This  possible  reduction 
in  size  and  weight  can  be  exploited  in  the  propulsion  of  hydrofoil 
craft.   Present  hydrofoil  propulsion  schemes  all  have  serious  salient 
disadvantages  that  have  stimulated  the  search  for  a  better  propulsion 
scheme  for  hydrofoil  craft.   This  thesis  proposes  the  use  of  a  novel 
dual  armature  superconducting  induction  motor.   The  motor  would  be 
encased  in  a  water  tight  pod  and  placed  on  the  end  of  a  foil  strut 
and  directly  drive  a  super-cavitating  propeller.  A  mathematical 
model  is  developed  and  the  size  and  parameters  of  a  21000  horsepower 
motor  are  determined.  The  resulting  model  is  statically  and 
dynamically  analyzed. 
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I  INTRODUCTION 

In  conventional  displacement  type  ship  hull  forms,  the  power  avail- 
able within  the  hull  can  "be  readily  converted  to  propulsive  thrust  by 
conventional  direct  drive  propellers.  However,  in  a  hydrofoil  type  ship 
hull  form,  the  mechanism  of  converting  the  available  power  within  the  hull 
to  propulsive  thrust  is  no  longer  a  simple  matter.  When  the  hydrofoil 
vehicle  'flies'  (becomes  foil  borne),  the  hull  is  lifted  out  of  the  water 
and  conventional  types  of  propulsion  are  no  longer  feasible.  A  simple 
direct  drive  straight  rotating  shaft  with  a  propeller  on  the  end  will 
not  suffice. 

An  assortment  of  alternative  hydrofoil  propulsion  schemes  have 
been  proposed  or  developed* 

1.  The  water  jet  system  scoops  up  water  at  an  inlet  in  a  forward 
foil  strut  and  ejects  this  water  at  high  velocity  out  the  stern, 
thus  providing  the  propulsive  thrust.  This  arrangement  requires 
large  amount  of  internal  hull  volume  for  water  ducts  and  pump 
machinery  and  places  large  amounts  of  water  weight  inside  the  hull. 
Additionally,  this  method  has  a  low  efficiency  at  low  velocities. 

2.  Another  scheme  utilizes  a  mechanical  gear  train  containing 
right  angle  gears  to  drive  propellers  mounted  on  the  foil  struts. 
This  scheme  suffers  from  a  very  high  initial  cost  and  a  low  mean 
time  between  failure,  hence  unreliability  and  high  maintenance  costs. 

3.  Air  propellers,  turbo-props,  and  other  aerodynamic  propulsors 
all  have  los  efficiencies  in  the  speed  ranges  anticipated  for 
hydrofoils  in  the  near  future. 
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k.       Electric  motors  can  be  placed  in  pods  on  the  end  of  the  foil 

struts  and  drive  propellers  directly.  This  method  eliminates  the 
objections  of  the  preceding  schemes,  but  introduces  a  new 
objection  of  its  own.  The  pod  of  the  size  required  for  a  conventional 
motor  would  produce  extremely  high  levels  of  hydrodynami c  drag  in 
the  speed  ranges  of  the  hydrofoil  vehicles. 

While  all  of  the  above  methods  of  propulsion  are  feasible  and  can  be 
utilized,  there  is  much  needed  improvements  before  hydrofoil  vehicles 
can  compete  economically  with  the  displacement  type  vessels.  Much 
research  and  development  is  being  done  to  improve  on  each  type  of  propul- 
sion system.   In  particular,  the  motor  in  the  submerged  pod  system  has 
been  vigorously  pursued  with  emphasis  on  reducing  the  required  size  of  the 

motor  through  the  use  of  super  conduction  technology.   '  '  '  ' 

The  first  attempt  at  applying  super  conducting  materials  to  motors 

(7) 

utilized  Faraday's  homopolar  concept.     It  has  the  characteristics  of 

low  DC  voltages  with  high  current  densities.  *^'     The  high  DC  current 
densities  of  the  homopolar  devices  complicates  the  power  transmission  and 
switching  problems  and  requires  some  sort  of  liquid  brushes  in  the  motor. 
Several  different  configurations  using  the  homopolar  concept  are 
discussed  in  reference  (4). 

The  recent  trend  in  super  conducting  machines  has  been  toward  the 
use  of  AG  synchronous  devices.  While  AG  machines  lack  the  high  starting 
torques  which  are  characteristic  of  DC  machines,  they  have  a  superior 
reliability  in  a  marine  environment  and  the  power  transmission  and 
switching  problems  are  much  less  than  for  DC  machines.   The  AC  super 
conducting  machines  most  considered  are  of  the  type  with  a  stationary  three 
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phase  armature  (stator)  and  a  revolving  DC  field  (rotor)  consisting  of 

super-conducting  winding  in  a  Dewar  bottle.^  *->*'' 

A  recent  development  in  AC  super-conducting  machines  is  a  family  of 
machines  which  utilizes  two  armatures  and  a  super-conducting  field 
winding.   '  The  field  winding  and  one  armature,  which  may  or  may  not  be 
wound,  rotate,  and  the  other  armature  is  fixed  (stator).  One  motor  config- 
uration utilizes  a  superconducting  field  which  is  free  to  rotate  and  is 
mechanically  unconstrained  except  for  its  own  inertia,  it  responds  only 
to  the  applied  electrical  torques  and  its  inertia.  Both  armature  windings 
are  wound  in  a  three  phase  manner,  one  armature  is  fixed  (stator)  and 
receives  the  input  electrical  power,  and  the  other,  positioned  between  the 
stationary  armature  and  the  rotating  field,  develops  the  mechanical  output 
torque  in  the  same  manner  as  a  wound  rotor  induction  motor  does.  This 
latter  configuration  is  the  topic  of  this  thesis.  It  will  be  mathematically 
developed  and  applied  to  the  propulsion  of  a  hydrofoil  vehicle.  Its  static 
and  dynamic  characteristics  will  be  presented.  The  resulting  size  of  the 
motor  will  demonstrate  that  it  is  compact  and  will  significantly  reduce  the 
hydrodynamic  drag  of  the  pod  to  an  acceptable  level. 


s 
II   PHYSICAL  DESCRIPTION  OF  THE  DUAL  ARMATURE  MOTOR 

The  motor  to  "be  studied  by  this  thesis  utilizes  a  four  pole 
rotating  field  that  is  mechanically  unconstrained,  except  for  its  own 
moment  of  inertia.  This  field  will  be  of  the  super  conducting  type  and 
will  be  encased  in  a  Dewar  bottle  and  cooled  by  liquid  helium.  It  will 
be  excited  by  a  DC  current  source.  In  general,  the  design  of  the  field 
will  be  similar  to  the  field  winding  developed  at  MIT.  ■"  This  field 
will  be  capable  of  producing  conventionally  unachievable  magnetic  flux 
densities  due  to  the  large  current  densities  that  are  made  possible  by 
the  super  conducting  material.  The  superconducting  winding  will  hence 
forth  be  referred  to  as  the  field. 

Surrounding  the  field  and  attached  to  it  is  a  damper  shell.  This 
shell  will  provide  the  start  up  torques  and  damn  out  dynamic  oscillations. 
The  damper  shell  and  the  field  together  will  hence  forth  be  referred  to  as 
the  rotor. 

Surrounding  the  rotor  will  be  a  three  phase,  four  pole  winding. 
This  winding  will  contain  no  magnetic  core  material  and  will  contain 
only  enough  steel  or  other  material  to  provide  mechanical  support  and 
rigidity.  It  will  be  free  to  rotate  relative  to  the  rotor  and  to  the 
motor  frame.   It  will  be  attached  to  a  shaft  from  which  mechanical  power 
may  be  extracted  and  supplied  to  the  load,  in  this  case  a  propeller.   The 
electrical  terminals  of  the  winding  will  be  brought  outside  the  motor 
through  slip  rings.   This  winding  will  hence  forth  be  referred  to  as 
the  shell. 

The  rotor  and  shell  are  surrounded  by  yet  another  three  phase  four 
pole  winding  that  contains  no  magnetic  core  material.   It  is 
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rigidly  attached  to  the  motor  frame  and  contains  only  enough  steel  to 

provide  mechanical  support  and  rigidity.  This  winding  will  hence  forth 
be  referred  to  as  the  stator. 

Surrounding  the  stator  is  a  shield  made  of  a  magnetic  core  material. 
Its  purpose  is  to  con-tain  the  magnetic  fluxes  within  the  motor  and 
provide  a  more  efficient  return  path  for  the  fluxes.  If  the  shield  were 
not  present  the  rotating  flux  field  would  extend  "beyond  the  motor  and 
possibly  induce  hazardous  voltages  in  the  surrounding  areas  which  could 
present  a  hazard  both  to  equipment  being  operated  in  the  vicinity  and 
to  man. 

The  stator  is  now  surrounded  by  a  water  tight  enclosure,  called  the 
pod,  that  will  keep  the  motor  dry  and  minimize  the  hydrodynamic  drag. 
The  magnetic  shield  could  be  a  part  of  the  pod  body. 

In  operation,  three  phase  electrical  power  will  be  provided  to 
the  stator  winding  which  will  produce  a  revolving  magnetic  field.  The 
rotor  will  lock  into  synchronism  with  the  stator  field  and  in  effect 
become  a.  virtual  magnetic  core  by  producing  a  rotating  magnetic  field 
far  greater  than  could  the  air  core  stator  by  itself.  The  shell  will 
act  as  a  wound  rotor  induction  armature  under  the  influence  of  the 
rotating  magnetic  field  and  drive  a  super  cavitating  propeller  directly 
as  shown  in  figure  (1).  Figure  (2)  shows  the  geometric  configuration 
of  the  field,  damper,  shell  stator  and  shield. 
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III  MATHEMATICAL  MODEL  OF  THE  MOTOR 

The  motor  is  modeled  electrically  as  shown  schematically  by  figure 
(3).  Figure  (3)  shows  the  axes  of  the  various  windings.  The  rotor  has 
a  direct  axis  (in  the  direction  of  positive  flux  lines  of  the  field)  and 
a  quadrature  axis  which  leads  the  direct  axis  by  ninety  degrees.  The 
direct  axis  is  at  an  angle  9  with  respect  to  the  stator  'A*  phase 
axis  as  shown  in  figure  (3)» 

The  shell  'A*  phase  axis  is  at  an  angle  jzf  with  respect  to  the  stator 
•A*  phase  axis. 

The  various  mutual  inductances  are  assumed  to  be  sinusoids  dependent 
upon  the  angular  orientation  of  the  axes  involved.  The  mutual  inductances 
are  maximum  when  the  axes  are  aligned,  zero  when  they  are  ninety  degrees 
apart,  and  negative  (minimum)  when  they  are  one  hundred -eighty  degrees 
apart.  Equation  set  (1)  gives  the  flux  linkages  for  the  various  circuits 
when  the  currents  in  all  circuits  are  considered.  Equation  set  (2)  gives 
the  voltages  at  the  terminals  of  the  various  circuits  when  all  currents  in 
all  circuits  are  accounted  for  and  the  circuit  relation  V.  ■=  pX.  +  R.  I. 
is  utilized.  The  symbol  p  is  a  operator  representing  the  time  derivative 
d/dt.  All  other  symbols  utilized  are  idetified  in  the  Glossary  of  Terms, 
Table  (1). 

The  rather  formidable  equation  sets  (1)  and  (2)  can  be  significantly 
simplified  by  performing  a  Park's  transformation  on  them.   The  Park's 
transformation  equations  are  given  here  as  equation  set  (3)»  When  the 
transformation  of  equation  set  (3)  is  applied  to  equation  sets  (1)  and 
(2),  equation  sets  (5)  and  (6)  result. 
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Table  I    Glossa.ry  of  Terms 

Note:  where  a  dual  assignment  of  usage  is  indicated,  the  context  should 

make  it  clear  to  the  reader  the  intended  usage  of  the  symbols. 

Symbols 

I     Current 

L     Self  Inductance 

1     length 

M  Mutual  Inductance 

N     Coil  Turns/  angle  as  defined  in  figure  (  3  ) 

P     Permeance 

p  Differential  Operator  d/dt 

R     Resistance/  Radius 

s     Slip 

T     Torque 

t     time 

U     Dummy  Variable  used  in  the  Park's  Transformation  Equation 

V     Voltage 

W     Angular  Velocity 

X     Reactance/  Ratio  of  Reactances 

9     Angle  Defined  in  Figure  (  3  ) 

X     Flux 

TT       3.TT^.... 

<j>  Angle  Defined  in  Figure  (  3  ) 

^  Per-Unitized  Flux 
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Table  I    continued 

Subscripts 

a  Phase  a/  Stator  Armature 

b  Phase  b/  Base  Quantity 

c  Phase  c 

d  Direct  Axis 

f  Field 

i  Arbitrary  Index 

j  Defined  in  Text 

k  Defined  in  Text 

1  Defined  in  Text 

m  Defined  in  Text 

o  Zero  Sequence/  Synchronous  Frequency/  Defined  in  Text 

p  Per-Unitized  Quantity 

r  Rotor 

s  Shell 


0.0. 
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Figure  (3) 
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Equation  set  (4)  can  now  be  obtained  directly  from  equation  sets  (5) 

and  (6). 

A  further  simplification  may  be  made  to  equation  set  (5)  by 
per  unitizing  it  in  accordance  with  reference  (3).  First,  equation  set  (5) 
is  re-written  as  equation  sets  (7),  (8)  and  (9). 

Then,  the  ordinary  variables  in  equation  sets  (7),  (8)  and  (9)  are 
normalized  in  accordance  with  reference  (3)  by  dividing  each  variable 
quantity  by  a  corresponding  base  quantity.  The  results  are  shown  as 
equation  sets  (10),  (11)  and  (12). 

It  is  highly  desirable  that  the  resulting  per-unitized  mutual 
inductances  exhibit  the  property  of  reciprocity,  that  is  M.  .  =  M. . . 
This  imposes  the  following  restrictions  on  the  base  quantity  selection. 

2.  v   1   =  2.  y   1   =  v   1   =v   1 

2  ab  ab   2  sb  sb    rb  rb    fb  fb 

Additionally,  the  following  conditions  on  the  base  quantities  are 
arbitrarily  imposed  to  simplify  the  per-unitized  equations. 

„M   W  I  ,    .,  M   WI,    M  _  W  I  „    (  L  -  M   )  W  I  , 

„  _  3  as  o  sb  _  3  ar  o  rb  _  af  o  fb  =  N   a    aa  /  o  ab 

a  2."        V  .  V  .         V  ,  V  , 

ab  ab  ab  ab 

Thus,  at  this  point,  six  constraints  have  been  placed  on  the  eight  base 

quantities.  This  leaves  only  two  base  quantity  parameters  free  to  be 

chosen  arbitrarily.  Once  these  two  more  conditions  are  specified  the  base 

quantities  will  be  completely  determined.  Further  consideration  of  base 

quantities  is  given  in  appendix  (H). 

Equation  sets  flo),  (11)  and  (1?.)  may  now  be  written  as  equations  sets 

(13),  (1*0  and  (15). 
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Appendix  (a)  shows  that  X  =  X   «=  X  _  and  X  =»  X  _  so  equation 

S      S3?      SI        3T      2?I 

sets  (i3)i  (14)  and  (15)  can  be  re-written  as  equation  sets  (16),  (17) 
and  (18)  respectively. 

If  the  base  quantities  utilized  to  obtain  equation  sets  (10),  (11) 
and  (12)  are  applied  to  equation  set  (4),  equation  set  (19)  results, 

Nov; j  if  equation  set  (19)  is  solved  for  the  time  derivatives  of 
the  fluxes,  equation  set  (20)  results. 

Equation  sets  (16),  (17)  and  (18)  can  be  inverted  so  that  fluxes 
are  the  dependent  variables  and  currents  are  the  independent  variables. 
The  resulting  equations  are  shown  as  equation  sets  (21),  (22)  and  (23). 
where,  X  (  X  -  X  ) 

Y   =   S      S      r 

j   X  (  X  -  X  j 
J    a  x  a    s  ' 


\m\ 

'  X     -  X     ) 

k     s         r 

;  x   -  x    ) 

'a         s 

xl  =  - 

;  x  -  x   ) 

'a         r 

;  X  -  X  ) 
k     a         s 

X    -  ■ 
m 

[  X  -  X  ) 
k      s         r 

;  xr  -  Xf  ) 

w 


Ao     X  -X  ) 

N  r    s  ' 

Equation  sets  (21),  (22),  (23)  and  (20)  can  be  combined  to  obtain 
equation  set  (24). 

The  direct  and  quadrature  axis  voltages  of  the  stator  armature  are 
given  by: 


V  .  =  V  COS  (  6  )   and  V 
adp     a     N  aqp 


-  V  SIN  (  0  ) 
a     v   ' 


adp       W     "    ad  W  aq         ap     adp 

V         -  J2_¥        +  Lj£l^        +  R       ! 
aqp       W         aq  W  ad         ap     aqp 


V         -  $-<£>        +  R       I 
aop       W         ao         ap     aop 

0         -  %-W    ,   -  i-J—l'i'  +  R       I   . 
W         sd          W            sq         sp     sdp 

0         -  £-¥         +  ^-W     .  +  R       I 

W         sq  W  sd         sp     sq 


0         =  £-¥         +  R        I 

W         so         sp     sop 
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£- *»     ,  +  R       I   . 

W         rd         rp    rdp 


0         =  $_tj>        +  r       i 

WQ       rq         rp     rqp 


Vf  p     3  Wq  *  f     +  Rf  p  If  p 


Equation  Set  (  19  ) 
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*      ad         o     adp       K*  '       aq         o     ap     adp 
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r      f  o     fp        o     fp    fp 


Equation  Set  (  20  ) 
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ad         o     adp       v  -^    '       aq         ap     o  x     k     sd         j       ad' 


piji        =WV         +  (  p6  )  *&     .  +  R       X     (  X,  4?  X .  »£»        ) 

r       aq         o     aqp       v  -^    '       ad         apovksq-j.aq' 

W     R 

p¥      -  w    v       -    °    a£*? 

r       ao         o     aop         X  ao 

r  ao 

p  ^         =  (  pN  )  *         +  R       X(X.^,-X..'*'+¥,) 
■^       sd       v  -^    '       sq         sp     o  v    k       ad         lsd  rd  ' 


p  *  =  -(pN)«S.+R       X     (  X,   ^         -  X-  «*         +  *£»        ) 

^       sq  v  ^    '       sd         sp     o  x     k       aq         1       sq  rq  ' 

p*i»  =  R       X    f*?     ,    -  (  1  +X    )^,+X    *E»  1 
-^      rd         rpoLsdv  m  '       rd         mfj 


^       rq         rp     o  v      sq  rq  ' 


W     R 

p*    -  w  v     -  °  ap>p 

•^       so         o     sop        X  so 

so 


p  *?_     -  W     V„  +  R^     X     X     ( »F     ,    -  SE»     ) 
^      f  of         fpom^rd  f  ' 


Equation  Set  (  2k  ) 
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The  dynamic  mechanical  equations  of  the  shell  and  rotor  in  per 

unit  form  are:     p2(  0  )  -  W  T  /z   H 

o  rp    r 

and      p2(6-N)=  W  (T  -T,  )/Z   H 
*   v   '   ov  sp  lp7/   s 

Appendix  (b)  provides  the  following  expressions: 

X 

T   =  77^  (  »5»    *i»    -  *•  ,  *&         ) 

rp   W  v   rd    sq     sd    rq  ' 

X 

T   =  —  (  «5*    vjj    -  ^    HP    )-T 

sp   W  v   sd    aq      ad    sq  '    rp 

These  equations  along  with  the  equations  of  equation  set  (24)  form 

a  complete  set  of  differential  equations  that  mathematically  describe^ 

the  motor's  model.  It  is  assumed  in  this  thesis  that  all  loadings  and 

transient  states  will  "be  symetrical  and  balanced.  Thus,  *f   and 

ao 

*F   will  he  assumed  to  be  zero.  Additionally,  *■£■»  _  will  be  considered 
so  r 

a  constant.   This  being  the  case,  the  given  differential  equations 
reduce  to  a  ninth  order  non-linear  set  of  differential  equations. 
These  equations  can  be  solved  once  all  the  constant  parameters  are 
known.  This  solution  is  rather  straight  forward  for  the  steady  state 
condition,  although  a  computer  minimizes  the  computational  effort. 
The  transient  state  defies  a  hand  solution.   The  digital  computer 
can  be  used  to  simulate  the  transient  state  once  a  set  of  initial 
conditions  has  been  determined. 

References  (8),  (9),  (10)  and  (12)  provided  valuable  background 
material  that  aided  in  the  development  of  the  equation  in  this  section. 
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IV  Determination  of  Machine  Parameters 


The  size  of  the  motor  is  determined  in  Appendix  (  C  )  and  the  circuit 
parameters  are  determined  in  Appendix  (  D ).  The  results  are  summarized  in 
Table  (  2 ) . 


Table  2    Machine  Size  and  Circuit  Parameters 
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Inner  Radius  of  Field 
Outer  Radius  of  Field 
Average  Radius  of  Damper 
Inner  Radius  of  Shell 
Outer  Radius  of  Shell 
Inner  Radius  of  Stator 
Outer  Radius  of  Stator 
Inner  Radius  of  Shield 
Outer  Radius  of  Shield 
Electrical  Length 

R 

ap 

R 
sp 


rp 


R 


fP 


3 


4.0  inches 

6.5   inches 

7.0   inches 

7>5   inches 

10.5  inches 

11.0  inches 

13*5  inches 

14.5  inches 

20.0  inches 

7.4  feet 

0.20 

0.588 

26.87 

0.0 

7.799 

3.685 

4.685 

10.08 

0.738  seconds 
0.0738  seconds 


(super-conducting  winding) 
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V  Steady  State  Characteristics 

When  the  motor  is  running  in  steady  state,  the  governing 
differential  equations  of  equation  set  (  2k  )  can  be  reduced  by  setting 
all  time  derivatives  of  the  fluxes  to  zero.  Additionally,  the  torque 

on  the  rotor  must  be  zero  and  this  directly  implies  <P   =  SP   ■  0  . 

J        r  rq     sq 


Slip,  s,  is  defined  defined  in  the  conventional  manner;  that  is 


W  -  W 

s  =  - where  W  =  p  (  9  )  »  W   and  W  =  p  (  <f>  ) 

r 

WQ  -  p  (  fi  ) 

in  the  steady  state,  then  s  ■ .  As  can  be  seen  from 

o 

figure  (3)e6,=  N  +  ^  sop(N)=W  -  p  (  jl  )  ,  so  finally  p  ( N)  »  s  W  . 

Taking  equation  set  (  24  )  and  setting  the  time  derivatives  of 

fluxes  to  zero,  letting  *5*   =  "P   =  0  ,  and  using  the  relations 

pN  =  s  W  ,  p9  =  W  and  I„  =  V „  /  R„  ,  equation  set  (  25  )  results. 

The  torque  equation  of  the  shell  reduces  to, 
X 
sp    W     sd   aq  and  by  using  the  fourth  equation  of 

equation  set  (  25  )»  the  expression  for  torque  can  be  written  as, 

R  X    q 

sp         s     v       W 


R  X    * 

T       ■  _§£  (   _2 a£_  n2 

o 


W  /  s  T 

Solving  for    ^P        yields        *P        =  =2.      /  s? 

*  aq  J  aq       R      ,  /        R 

^*  ^         spy  sp 

Assuming  that  slip,   s,   stator  torque,   T     ,   and  field  current,   I„, 

sp  I 

are  the  independent  variables,  then  equation  set  (  25  )  can  be  solved 
for  the  fluxes  and  voltages  which  result  in  equation  set  (  26  ). 

The  digital  computer  is  a  valuable  tool  in  evaluating  the  expressions 
of  equations  set  (  26  )  numerically.  Equation  set  (  26  )  and  the  machine 
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0  =  X,  «P  .  -  X,  *  ,  +  ¥  , 

k   ad    1  sd    rd 


R   X 

0  -  s  «P  .  +  -~^— -  X.  ¥ 

sd     W     k   aq 
o       ^ 


sd   v     m  '     rd    m   f 


o-i,  +  5ti(»      *,) 


Equation  Set  (  25  ) 


35 


W  /  s  T 

^aq       X      ./       R 

^         ON/  SP 

R       X  W     I- 

ad         s  W       x  1'       aq       X.    X 

o  ^         k     o 


qa        =0 
sq 

R       X 
vj,        -  -  _^_o     x    ^ 

sd             s  W           k       aq 
o                  ^L 

<P        =0 
rq 

rd           sd 

X 
o 

W 

*  f             rd 

+       °       I 

X     X     xf 
o     m 

V    „      =   _  cu 

R       X 

.      ap     o  / 

adp  aq  W         v      j       ad         k       sd 


R       X 

V         ■**     .  +     5^     °  x-  * 

aqp  ad  W  J       aq 


v+  =x/(V  A    f  +  (V        )2 
t      v       adp'  a<lP 

Equation  Set  (  26  ) 
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parameters  listed  in  Table  (  2  )  along  with  the  load  characteristics 

determined  in  appendix  (  E  )  were  injected  into  an  IBM  model  370 

computer  use  the  WATFIV  program  listed  in  appendix  ( G  ). 

The  results  of  the  computer  analysis  are  shown  in  the  results 

section  as  figure  (  4  )  through  (  15  )• 


3? 

VI  Transient  State  Characteristics 

Equation  set  (  Zh  ),   table  (  2  )  and  the  shell  and  rotor  torque 
expressions  of  Appendix  (  B  )  provide  the  mathematical  discription 
of  the  transient  state.  Given  these  equations  and  a  complete  set  of 
initial  conditions  and  driving  functions,  the  system  is  deterministic. 
The  initial  conditions  are  determined  in  Appendix  (  F  ). 

The  transient  state  was  simulated  in  four  distinct  evolutions: 

1 .  Rotor  start  up  with  R   =  infinity 

2.  Shell  start  up  with  rotor  at  full  speed  and  at  various  values  of 

R 
sp 

3.  Crash  back. 

h.  Three  phase  fault  at  the  stator  armature  terminals. 

The  simulations  were  all  conducted  using  an  IBM  model  370  digital 
computer  and  the  IBM  supplied  program  CSMP  (Continuous  Systems  Modeling 
Program).  The  program  listings  are  shown  in  Appendix  (  G  ). 

The  rotor  start  up  simulation  assumed  that  the  shell  circuit  is 
open  and  in  effect  not  there.   This  was  equivalent  to  starting  a 
conventional  single  armature  machine.   The  results  of  this  start  up 
ar  shown  in  figure  (  16  ). 

The  shell  start  up  simulation  was  conducted  by  first  having  the 
rotor  running  in  steady  state  then  placing  various  values  of  finite 
resistance  in  the  shell  circuit.   The  results  of  this  simulation  are 
shovrn  in  figure  (  17  ). 

The  crash  back  simulation  first  assumed  a  steady  state  full  power 
condition  for  the  motor,  then  the  phase  sequence  of  the  stator  terminal 
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voltage  were  reversed.  The  results  of  the  crash  back  simulation  are 

shown  in  figure  (  18 ). 

The  three  phase  fault  simulation  again  assumes  a  steady  state  full 

power  condition  for  the  rotor,  then  the  three  terminal  leads  of  the  stator 

armature  are  electrically  shorted.  In  effect  V   =  V  ,  =  V   =0. 

aa    ab    ac 

The  results  of  the  three  phase  fault  are  shown  in  figure  (  19  )• 
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VII  Results 

Figure  (  4 )  through  figure  (  15  )  show  the  steady  state  analysis 
results.  In  these  figures,  s,  I   and  T   are  the  independent  variables. 

Figures  (  4  )e  (  5  )  and  (  6  )  show  the  terminal  voltages. 

Figures  (  7  )»  (  8 )  and  (  9  )  show  the  terminal  currents. 

Figures  (10),  (11)  and  (12)  show  the  power  factor. 

Figures  (13)»  (14)  and  (15)  show  the  rotor  angle. 
Figure  (  16  )  through  figure  (  19  )  show  the  angular  velocities  of  the 
shell  in  the  transient  state. 

Figure  (  16  )  shows  the  rotor  start  up  transient. 

Figure  (  17  )  shows  the  stator  start  up  transient. 

Figure  (  18  )  shows  the  crash  back  transient. 

Figure  (  19  )  shows  the  three  phase  fault  transient. 
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VIII   DISCUSSION  OF  RESULTS  AND  CONCLUSIONS 

Figures  (4)  through  (15)  demonstrate  that  the  dual  armature 

motor  can  be  characterized  by  a  set  of  steady  state  condition  curves. 

These  curves  can  assist  in  the  rational  selection  of  an  operating 

condition  for  the  motor.   The  data  supplied  by  these  curves  can  be 

cut  many  different  ways;  for  instance,  I„  or  T   could  be  considered 

fp     sp 

as  dependent  variables  and  V,  or  I,  could  be  treated  as  independent 
variables.   In  cutting  the  data  different  ways  additional  data  points 
may  be  desirable,  in  this  case  the  WATFIV  program  in  Appendix  (G)  could 
be  used  if  the  appropriate  parameters  were  adjusted. 

Figures  (16)  through  (19)  provide  the  transient  state  results. 
Figure  (16)  shows  that  the  rotor  on  start  up  under  goes  large  angular 
velocity  excursions  with  the  attendant  large  acceleration  forces.  At 
one  point  the  rotor  actually  runs  bqr?T<wards.  Eventually  the  rotor 
settles  down  to  its  mechanical  synchronous  speed  of  1800  RPM. 

Figure  (17)  shows  that  for  the  shell  start  up  there  is  some  critical 
parameter  R   between  10  pu  and  20  pu  where  the  shell  starts  up  in  a 

stable  fashion  at  the  fastest  possible  rate.  For  values  of  R   less  than 

r  sp 

this  critical  value,  the  rotor  tumbles  or  pulls  out  of  synchronism  with  the 

stator  field  and  results  in  a  loss  of  starting  torque  and  the  introduction 

of  a  pulsating  torque.  For  R   =  O.588  pu  (its  minimum  value)  the  motor 

sp 

rotor  does  not  start  at  all  but  oscillates  about  the  zero  point. 

Figure  (18)  shows  the  crash  back  W   curve.  The  results  of  this 
curve  beyond  1,8  seconds  are  questionable.   The  computer  simulation  appears 
to  have  gone  astray  at  this  point,  perhaps  due  to  too  large  of  a  step 
size  in  the  simulation  integration  routine  or  due  to  the  proximity  of  W 
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to  zerd  The  rotor  and  shell  both  move  at  approximately  the  same  velocity 
in  this  maneuver. 

Figure  (19)  shows  the  three  phase  fault  at  the  stator  armature 
terminals.  Again  the  computer  simulation  appears  to  go  astray  at  about 
3.7  seconds  so  results  beyond  that  point  are  questionable. 

The  preliminary  results  of  this  thesis  indicate  that  the  dual 
armature  super  conducting  motor  is  indeed  a  candidate  for  hydrofoil 
propulsion  and  warrants  further  investigation. 
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IX   RECOMMENDATIONS 

Additional  areas  that  could  be  investigated  using  the  equations 
developed  in  this  thesis  nay  he: 

1.  Determine  steady  state  conditions  as  a  function  of  R   .  The 

sp 

present  investigation  considered  R   =  O.588  pu  in  the  steady  state 
condition. 

2.  Write  a  computer  program  that  will  accept  machine  dimensions  as 
input  and  return  machine  parameters  for  use  in  the  steady  state 

or  transient  analysis. 

3.  Instead  of  holding  SF  _  constant  as  done  in  this  thesis  in  the 
transient  state  investigation,  hold  I„  constant. 

h-.       Investigate  current  densities  in  the  various  windings  during 
the  transient  periods  to  see  if  they  exceed  tolerable  levels.  Devise 
a  scheme  for  controlling  excess  currents  during  transients. 

5.  Determine  the  static  stability  limit  of  rotor  angle. 

6.  Refine  load  characteristics. 

7.  Vary  machine  dimensions  to  optimize  the  design.  Criteria  for 
optimizing  could  also  be  developed  as  a  side  issue. 

8.  Introdue  hydrofoil  dynamics  into  the  model. 

9.  Investigate  speed  control  by  varying  W  the  electrical  fre- 
quency of  the  prime  mover. 

Most  of  these  subject  areas  would  apply  at  an  under  graduate  project 
level. 
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Appendix  A    Shell  Constraints 

Equation  sets  (  13  )»  (  1^  )  and  (  15  )  contain  parameters  that  are 
not  independente  The  following  development  assumes  that  flux 
leakages  between  adjacent  circuits  is  negligible . 

Assuming  that  the  stator,  shell  and  field  axes  are  aligned  so  as 

to  eliminate  the  sinusoidal  dependency  of  the  mutual  inductances ,  the 

following  relations  can  be  written. 

(L-M       )  =|n2  (   P  +  6P      ) 
vs         ss'2sx  ss   ' 


\  -  Nr2  (   P  +  6  £r   ) 


M 
as 

B 

N 
a 

N     1 
s 

^  P  +  6  Pas  ) 

M 

ar 

= 

N 
a 

N     ( 
r 

'   P+6P      ) 
ar  ' 

M 
af 

= 

N 
a 

Nf  1 

:  p  +  6  paf ) 

M 
sr 

= 

N 
s 

N     < 
r 

:  p  +  6Psr) 

Msf 

= 

N 
s 

Nf  1 

'  P  +  6  Psf  ) 

Mrf 

=: 

N 
r 

Nf  ( 

:  p  +  6  prf ) 

The  *  factor  appearing  in  the  first  equation  is  characteristic  of 
three  phase  windings. 

Now,  if  the  various  leakage  permeances  may  be  considered  negligible 
and  are  set  to  zero,  the  above  equations  may  be  used  to  derive  the 
following  relations. 
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2  af  v  a   "aa  ;   2  sf   2  as 

|  M   (  L  -  M    )=Jm    '  J  M 
2  ar  v  a    aa  '   2sr   2as 

K  „  L  -  M  _  M 
af  r    rf  ar 

These  equations  may  now  be  per-unitised  as  follows s 

2.    af  o  ab  _  s    'ss'  o  st  _  2  sf  o  sb  2.  as  o  ab 

2     v^  V  .  2  "   V  _      2    V  . 

fb  sb  fb  sb 

Q  M   WI,(L-M)WI,     0  M   W  I  ,   „  M   V  I  . 

2.     ar  o  ab  v  s    ss   o  sb     2,  sr  o  sb  2   a,s  o  ab 

2     V  .  V  ,     "    =  2     V  .     2     V  . 

rb  sb  rb  sb 

M  „  W  I_   L  W  I  .    M  ,  W  I _   M   W  I  , 

af  o  fb   r  o  ro  _  rf  o  fb   ar  o  rb 

ab        rb         rb         ab 

Comparing  these  relations  with  those  in  equation  sets  (  10  )  and 

(  13  )  it  can  be  seen  that  the  above  relations  reduce  to; 

X  -  X  . 
s    sf, 

X  -  X 
s    sr  and 

X  =  X  _ 

r    rf. 
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Appendix  B     Derivation  of  Torque  Expressions 

The  torque  on  any  winding  in  per  unit  is  given  "by 

Tip  =  '  ***  iq  Iidp   ^  id  Iiqp  '• 
Applying  this  to  the  stator  and  shell  windings  and  substituting  the 
expressions  for  current  from  equation  Sets  (21)  and  (22)  and  canceling 
terms  where  possible  yields: 

X  X 

T   -  °  K  (  ^  .  »P   -    f    ¥  ,  )   and 
ap    W    v   ad    sq      aq    sd  ' 

X 

T   =  rr^  (  *  ,  ^   -  ¥     «f  ,  )  -  T 

sp   W  v   sd   rq      sq    rd  '    ap 


Also  equilibrium  of  forces  requires  T   +  T   +  T   =  0  . 

ap    sp    rp 

Solving  the  above  equations  for  the  torque  on  the  rotor  gives 

X 

T   -  ~  (  5*    f.  -  ¥  ,  ¥    )  . 

rp   W  v    sq    rd      sd    rq 
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Appendix  G   Motor  Sizing 

References  (11)  and  (14)  contain  equations  for  the  power  rating  of 
an  air  core,  three  phase  synchronous  motor.  They  can  be  simplified  and 
adapted  to  the  present  application  by  assuming  the  following? 

a.  number  of  pole  pairs  =  2 

b.  winding  angle  of  field  =  120° 

c.  winding  angle  of  shell  -  60° 

d.  winding  angle  of  stator  =  60° 

e.  all  winding  electrical  lengths=  1 

f.  the  following  parameters  are  defined  as 

X  =  Rai/Rao  •  y  =  Rfi/Rfo  »  and  z  "  Rsi/Rso  '  The  radius 
dimensions  are  shown  in  figure  (C-i). 

The  resulting  power  rating  expression  is, 

P ■  -  12  /T  10-1  Wo  Ja  Jf  c2  R^  (1  -  y*)(  i  ) 
where  C    -  |  (  -  In  x  +  ^  "  x  ■)  (  -2£  f)      f 


-2       g  v       ^  A  ^  v  - 

sh 

B4-  -XaSIN(Y)  -^/l   -xa  C0S2(Y)      , 


COS(y)  =  power  factor    , 

J-    G_  R£ 

Vo 


h-yfiZT      J;(li  ^  ( Re  }\ 


-s       ■%«—-■       4  +  g  (   g 


and         Co  «  1  (    x4  In  x  +  ^=^-     +  i    1    =  x     )2(   ^o   ^  ) 


sh 
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Figure  (  C  - 1  ) 
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The  current  densities  were  chosen  based  on  the  values  used  in 

reference  (11)  and  are;     J  =  2.50  x  10  amps /meter  t 

Q  p 

and     J„  «=  1.25  x  10  amps/meter  , 
Synchronous  electric  frequency  was  chosen  to  be  W  =  377  radians/sec. 
Power  factor  was  chosen  to  be  PF  =  0.8  * 

The  power  was  chosen  based  on  the  requirements  specified  in  reference 
(4)  and  is  P  =  2lt000  horespower  per  motor « 

The  radii  selected  are; 

R„.  =  4.0  inches  R_  =  6,5   inches 

fi  fo 

R  .  =  7,5   inches  R   =10.5  inches 

si     J  so     J 

R  .  =11.0  inches  R   =  13.5  inches 

ai  a-o 

R   =7.0  inches  Rgh  =  14.5  inches 

Finally,  if  the  above  assumptions  are  made  and  the  power  rating 
expression  is  solved  for  1,  the  result  is  1  -  7.4  feet. 
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Appendix  D   Circuit  and  Mechanical  Parameters  of  Motor 

References  (  11  )  and  (  14  )  contain  equations  for  determining  the 

inductances  (both  self  and  mutual)  for  an  air  core,  three  phase, 

synchronous  motor.  These  equations  are  infinite  series  summations j 

however,  the  lead  term,  n=l,  dominates  and  all  other  terms  are  small. 

These  equations  can  be  simplified  and  adapted  to  the  present  application 

by  assuming,  in  addition  to  the  assumptions  of  appendix  (c),  that? 

a.  the  leading  term  of  the  infinite  series  dominates, 

b.  the  shield  is  made  of  magnetic  iron  and 

c.  the  theory  of  superposition  holds. 

The  resulting  inductance  expressions  are  listed  in  Table  (  D  -  1  ) . 
The  dimensions  are  the  same  as  discussed  in  appendix  (c)  and  shown  in 
figure  (  C-l  ).  If  these  dimensions  are  now  used  to  solve  the  expressions 
in  table  (  D-i  ),  table  (  D-2  )  results. 

If  the  expressions  of  table  (  D  -  2  )  are  per  unitized  using  the 
base  quantities  of  appendix  (  H  ),  table  (  D-3  )  results. 

The  definitions  of  X.,  X^  X.^  Xm  and  Xq  given  on  page  25  can  now  be 
evaluated  and  are;      X.  =  7.799 


Xk  =3.685 

X2  -  4.685 

X  -  10.08 
o 

Reference  (14)  provides  an  equation  for  R  and  is 

17  1  N  2 
R 


r   4  a    Rx 
s  t 


where  C^   =  o   tj,  t.  being  the  damper  shell  thickness  and  O   being  the 


Table  D-l       Inductances 
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Table  D~2       Inductances 


L     -   5.104  x  10"7  N  2  1 
a  a 

L     =  3*532  x  10"7  N  2  1 
s  s 

L     =  2.601  x  10"7  N  2  1 
r  r 

Lf  =  2.101  x  10°7  Nf2  1 


Table  D»3       Pe£  Unitized  Inductances 


X     =  9.089 

x^  =  19.23 
x   -  56.62 

r       -^ 

X_  -  112.4 
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conductivity  of  copper.  Assuming  O   =  3,6  x   10      ^meter  and 

t  =  0.11  inches,  R  =  k.k2   x  10  J   N   1.  Per  unitizing  using  the 

base  quantities  of  appendix  (  H  )  gives,  R   =  26.87  pu. 

In  general,  the  resistance  of  a  wound  three  phase  winding  is 

2 
given  by  4  N.   1. 

R,  =  ' 


1    (R  .2  -  R..2)  6   .  K  .  G. 

x  oi     li  '  wei  si  l 


assuming  the  conductivity  is  CJ  =  6  x  10  Ehos/meter  and  the  space  factor 

K  =  0.27  for  both  the  shell  and  the  stator  armature,  the  per  unitized 

resistances  are;     R   =  0.200  pu  and   R   =  0.588  . 

ap  sp 

The  resistance  of  the  field  winding  is  of  course  zero  because 
it  is  made  of  a  superconductor. 

An  initial  estimate  of  the  moments  of  inertia  for  the  rotor  and 

for  the  shell  with  shaft, propeller  and  entrained  water  included  were; 

J  =  120  slug  ft2    and  J  =  12  slug  ft  2 
s        -to  r        & 

For  H  defined  in  the  conventional  manner,  that  is, 

2 

J.  W, 

H  -   i  b 


2  P 
b 


H  =  0.738  sec  and      H  •=  0.0738  sec. 
Note:  the  computer  listings  of  appendix  (G)  show  H to  be  different 

from  the  above  values.  This  is  because  W  was  included  so  the  value  in 

o 

the  listings  is  377  times  smaller  than  the  actual  assumed  value. 
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Appendix  E  Load  Characteristics 


The  dynamics  of  the  motor  and  propeller  system  are 

considered  to  be  uncoupled  from  the  hydrofoil  dynamics  for  the 

time  span  of  the  simulations  since  the  time  constants  of  the  two 

systems  are  quite  different.   Thus  the  load  torque  characteristics 

are  a  function  of  W  alone, 

s 

For  the  shell  start  up, the  propeller  characteristics  are 


assumed  to  he 


T,  -  5  (w  /to  )  . 

lp   y  s   s'  o' 


For  the  crash  astern  and  the  three  phase  fault,  the  load 
characteristics  are  assumed  as  shown  below. 

2 


-r 


-3- 


-1-- 


sp 


/ 


/ 


sp 
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Appendix  F  Determination  of  Initial  Conditions 

An  operating  point  was  selected  for  the  full  power  steady  state 

from  figures  (k)   through  (15)  as  s=0.03l5,  I  =1.0  pu,  and  T  =1.0  pu. 

ip  sp 

Then  using  the  equations  of  equation  set  (26)  a  complete  set  of  initial 

condition  were  determined  and  are  listed  in  the  computer  listings  in 

appendix  (G). 

For  rotor  start  up,  it  was  assumed  that  T  =0.0,  W  =0.0,  W  =0.0, 
n  sp    *   r    f  s    ' 

and  9=0.0  .   The  remaining  initial  conditions  were  calculated  using 
Equation  Set  (26)  and  are  listed  in  the  computer  listings  in  appendix 
(G). 

Part  of  the  rotor  start  up  computer  simulation  routine  was  to 
print  out  the  final  values  of  all  the  time  variables  at  the  end  of 
the  simulation  period.  These  final  values  were  then  utilized  as  the 
initial  conditions  for  the  shell  start  up.  The  values  are  listed  in 
the  computer  listings  in  appendix  (G). 
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Appendix  G  Computer  Program  listings 
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Appendix  H  Base  Quantity  Selections 

As  mentioned  on  page  (  21  ),  two  remaining  base  parameters  are  free 
to  be  selected.  The  bases  selected  are; 

1)  ^  =  21 , 000  horsepower  and 

2)  Ifp  =  1.0  pu  when  Jf  =  1.25  x  108  amps/meter2. 
At  this  point  eight  conditions  have  been  specified  on  the  eight 

base  quantities  so  the  base  system  is  fully  constrained  and  all  base 
quantities  can  be  calculated.  They  are; 

Vab  =  l9'36  Na  volts  Iab  -  2.696  x  105  amps/Na 

lsb  =  ^M9   x  105  amps/N 
Irb  -  i.388  x  106  amps/Nr 
1^  =  2.217  x  106  amps/N 


Vsb 

= 

II.65 

N 
s 

volts 

Vrb 

= 

5.640 

N 
r 

volts 

Vfb 

= 

3.531 

Nf 

volts 
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